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I.  Introduction 

Time  Domain  Finite  Difference  Method  was  a  novelty  in  field  computations, 
because  it  had  the  potential  to  solve  many  hitherto  unsolvable  problems,  but  it  required 
so  much  memory  that  it  was  beyond  the  available  resources  for  many  research  labora¬ 
tories.  But,  that  limitation  is  no  longer  there.  Because  of  the  reasonable  availability  of 
the  Cray  machines.  It  is  therefore  essential  to  improve  the  time  domain  computational 
method  so  that  it  can  be  applied  to  practical  problems.  This  report  presents  the  accom¬ 
plishments  of  the  work,  supported  by  the  contract,  in  advancing  the  art  of  time  domain 
computation  both  in  theory  and  practice. 

II.  Time  Domain  Finite  Difference  (TDFD) 

* 

The  time  domain  finite  difference  method  as  presented  by  Yee  ,  uses  the  leap-frog 
method  to  solve  scattering  problems  is  not  boundary  conforming,  i.e.,  it  uses  steps  to 
approximate  curved  surfaces.  Besides  the  inaccuracy  incurred  as  a  result  of  the  step 
approximation,  the  method  is  not  satisfying  as  a  modeling  technique  of  fine  structures  for 
many  microwave  components.  It  is  also  rather  common  that  the  time  domain  technique 
of  Yee’s  is  applied  with  only  a  monochromatic  wave  as  excitation. 

A  scrutiny  of  many  of  the  prior  computations  has  led  us  to  believe  that  the  time 
domain  method  needs  some  honest  clarification  before  large  scale  applications  are  possi¬ 
ble.  One  issue  is  the  conforming  boundary,  and  the  other  is  the  absorbing  boundary  con¬ 
dition.  Both  problems  are  not  emphasized  by  those  who  promote  the  time  domain  tech¬ 
niques.  They  seem  to  be  irritations  but  not  road  blocks  to  the  applications  of  the  method. 
Yet,  without  putting  these  two  issues  to  rest,  it  is  almost  impossible  to  rely  on  the  time 
domain  method  for  dependable  results.  That  is  why  the  objective  of  this  project  is  to 
emphasize  the  development  of  the  time  domain  finite  element  method  and  the  local 
absorbing  boundary  conditions. 

III.  Time  Domain  Finite  Element  Method 

To  relax  the  restrictions  of  the  finite  difference  method,  we  use  an  interpolation 
technique  to  find  the  field  inside  an  element  from  the  surrounding  nodal  points.  The  dif¬ 
ferentiation  of  the  function  is  analytically  found  from  the  interpolation  function.  Since 
the  interpolation  is  not  limited  to  regular  mesh,  we  are  able  to  generalize  the  method  to 
conforming  meshes.  This  technique  is  fully  developed  during  the  course  of  this  project 
and  resulted  in  a  Ph.D  dissertation  and  several  conference  presentations  and  one  journal 
publication. 

IV.  Time  Domain  Absorbing  Boundary  Condition  at  Material  Interfaces 

The  absorbing  boundary  condition  is  a  necessary  pan  of  solving  radiation  problems 
in  any  differential  equation  approach.  However,  it  has  not  been  considered  a  serious 
problem  by  most  practitioners.  Ad  hoc  methods  have  been  used,  although  many  interest¬ 
ing  papers  have  been  published  in  that  subject  The  reason  is,  in  most  scattering  prob¬ 
lems,  those  ad  hoc  method  seem  to  give  reasonable  results. 

♦  _ 

Yee,  K.  S.  "Numerical  Solution  of  Initial  Boundary  Value  Problems  Involving  Maxwell’s  Equations  In 
Isotropic  Media,”  IEEE  Trans.  Ant.  and  Propagation,  Vol.  AP-14,  pp.  302-307, 1966. 


The  situation  is  quite  different,  however,  when  one  applies  the  time  domain  method 
to  scattering  by  buried  targets  or  to  modeling  of  microwave  components.  The  emphasis 
on  that  topic  in  this  project  has  also  resulted  in  one  Ph.D  dissertation  and  several  confer¬ 
ence  and  journal  papers. 

V.  Ph.D  Dissertations 

The  Ph.D  dissertations  supported  in  whole  or  in  part  by  this  project  are: 

1.  “Time  Domain  Computation  of  Electromagnetic  Wave  Scattering  by  the  Method  of 
Conforming  Boundary  Elements,”  by  Andreas  Costa  Cangellaris. 

2.  “Numerical  Modeling  of  Two-Dimensional  Time  Domain  Electromagnetic  Scatter¬ 
ing  by  Underground  Inhomogeneities,”  by  Chung  Chi  Lin. 

Abstract  of  the  above  dissertations  are  attached. 
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2.  “Point-Matched  Time  Domain  Finite  Element  Methods  for  Electromagnetic  Radia¬ 
tion  and  Scattering,”  IEEE  Trans,  on  Antennas  and  Propagation,  September  1987. 

Conference  Presentations: 

1.  “Time  Domain  Absorbing  Boundary  Condition  in  Lossy  Media,”  IEEE  AP-S/URSI 
Symposium  Digest,  pp.  7  -  10,  Boston,  MA,  June  1984. 
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ence  Digest,  p.  227,  Philadelphia,  June  1986. 
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ABSTRACT 

The  method  of  conforming  boundary  elements  is  developed  and  applied  for 
the  numerical  computation  of  transient  electromagnetic  wave  scattering  from 
composite  two  and  three-dimensional  structures  in  free  space  illuminated  by  an 
impulsive  wave. 

The  transient  scattering  problem  is  formulated  as  an  initial-boundary-value 
problem  governed  by  Maxwell’s  differential  equations  and  is  numerically  solved 
by  means  of  a  finite-element  time-domain  method.  Boundary  elements  are  intro¬ 
duced  to  conform  the  desirable  conducting  geometries  and  dielectric  interfaces. 
Additional  numerical  conditions  are  implemented  for  satisfying  the  radiation  con¬ 
dition  for  the  scattered  field  on  an  artificial  boundary  that  is  used  for  the  trunca¬ 
tion  of  the  finite-element  region.  The  convergence  and  the  stability  of  the 
numerical  method,  as  well  as  the  various  types  of  errors  introduced  by  the 
numerical  discretization,  are  investigated.  Stable  numerical  models  for  treating 
the  field  singularities  at  the  vicinity  of  edges  and  corners  are  also  discussed. 

Computations  are  carried  out  for  scattering  by  two-dimensional,  perfectly 
conducting  cylinders  of  arbitrary  cross-section  illuminated  by  a  plane-wave  Gaus¬ 
sian  pulse  at  normal  incidence.  The  results  are  in  good  agreement  with  those 
obtained  from  a  time-domain  integral  equation  approach.  Scattering  by  dielectric 
cylinders  is  also  computed  and  the  focusing  effect  of  the  dielectrics  is  illustrated. 
Finally,  the  method  is  applied  for  the  numerical  computation  of  scattering  by 
finite  conducting  cylindrical  structures  with  constant  cross-section.  The  structure 
can  be  more  complex  by  adding  conducting  fins  or  wires  to  it.  Available  experi¬ 
mental  results  provide  us  with  a  further  check  of  the  validity  of  the  method. 
Because  of  its  demonstrated  accuracy,  efficiency  and  versatility,  the  method  will 
be  a  powerful  tool  in  the  study  of  electromagnetic  scattering  and  penetration  of 
large  arbitrary  shaped  composite  targets. 


Numerical  Modeling  of  Two-Dimensional  Time  Domain 
Electromagnetic  Scattering  by  Underground 
Inhomogeneities 

h 

Chung- Chi  Lin 


ABSTRACT 

The  preseat  study  examines  the  feasibility  of  the  active  electromagnetic 
remote  sensing  of  elongated  subsurface  structures  such  as  hidden  archeological 
ruins,  tunnels,  conduits  and  underground  waterfiows.  The  theoretical  investiga¬ 
tion  of  the  scattering  from  a  buried  target  in  the  homogeneous  earth  is  performed 
in  the  presence  of  an  impulsive  incident  wave  excitation. 

The  initial-boundary-value  problem  formulation,  by  Maxwell’s  differential 
equations,  of  two-dimensional  scattering  b  numerically  solved  by  means  of  the 
point-matched  (collocation)  finite  element  method  in  the  time  domain.  The 
greatest  difficulty  encountered  in  this  modeling  process  b  in  satisfying  the  radia¬ 
tion  condition  of  the  scattered  field  at  the  boundary  of  the  finite  element  region. 
An  extension  of  the  free  space  absorbing  boundary  condition  allowed  to  account 
for  the  lossy  character  or  the  earth  and  the  presence  of  the  air-earth  interface, 
achieving  a  better  absorbtiou  of  the  outward  travelling  wave. 

Computations  an  carried  out  for  two  different  types  of  excitations:  a  plane 
wave  at  the  normal  incidence  and  a  line  source  above  the  ground.  The  time  his¬ 
tories  of  the  scattered  fields  along  the  ground  surface  are  presented  and  analysed 
through  a  parametric  study.  Those  results  demonstrate  the  definite  advantage  of 
the  fast  transient  sensing  method  over  the  use  of  the  more  conventional 
sinusoidal  excitation.  In  particular,  the  time  delay  of  the  specular  reflection  from 
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the  target  with  respect  to  the  ground  response  facilitates  the  identification  of  the 
inhomogeneity’s  presence. 
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Point-Matched  Time  Domain  Finite  Element 
Methods  for  Electromagnetic  Radiation 

and  Scattering 

ANDREAS  C.  CANGELLARIS,  member,  ieee,  CHUNG-CHI  LIN,  member,  ieee,  and  KENNETH  K.  MEI,  fellow,  ieee 


Abstract—  Direct  time  domain  computation  of  Maxwell's  diffeteatiai 
equations  will  soon  become  a  practical  technique  because  of  the 
availability  of  supercomputers.  The  principal  methods  nsed  in  time 
domain  computations  and  the  supporting  theories  are  presented.  The 
point-matched  finite  element  method  is  chosen  as  the  main  feature  of  this 
presentation,  which  indndes  the  discretisation  of  equations,  conforming 
mesh  genera  lion,  dielectric  and  metallic  interfaces,  anmericai  stabBty 
and  simulation  of  mdiatioo  conditions.  Namericni  results  of  scattering  of 
Gaussian  poises  are  presented  in  time  sequences. 


of  Electromagnetic  W«»*<  and  Applications,  VoL  1,  Jf».  3,  101-132,  1987 
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Unimoment  Method  for  Electromagnetic  Ware  Scattering 
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Abitract-A  simple  problem  in  2-dimensional  space  is  used  as  an  example  for  the  applica¬ 
tion  of  the  unimoment  method.  The  discos  lion  includes  fundamentals,  mesh  generation, 
loading  of  equations,  sparse  matrix  inversion  and  boundary  matching.  The  objective  of 
this  paper  is  to  help  uninitiated  readers  understand  and  apply  this  method. 


